Riemannian light cone from vanishing birefringence 
in premetric vacuum electrodynamics 
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We consider premetric electrodynamics with a local and linear constitutinen7vareve law for the 
vacuum. Within this framework, we find quartic Fresnel wave surfaces for the propagation of light. 
If we require (i) the Fresnel equation to have only real solutions and (ii) the vanishing of birefringence 
in vacuum, then a Riemannian light cone is implied. No proper Finslerian structure can occur. This 
is generalized to dynamical equations of any order. 

PACS numbers: 03.50.De, 04.80.-y, 41.20.-q, 03.30. +p 



I. INTRODUCTION AND MOTIVATION 

Recently, the physics of the electromagnetic field, with- 
out assumptions about the metrical structure of the un- 
derlying spacetime, has gained renewed interest. On the 
one hand, this general ansatz is needed for a proper inter- 
pretation of experiments testing Lorentz invariance. In 
such approaches it is not allowed to make assumptions 
about the underlying geometric structure, in particular, 
about a metric of spacetime. On the contrary, by using 
properties of the evolution of the electromagnetic field, 
one likes to establish the metrical structure of spacetime 
(here "metrical" may be more general than the ordinary 
Riemannian or Minkowskian metric) . The general struc- 
ture of Maxwell equations can serve as a test theory for 
searches for Lorentz violation in the photon sector 0. 
On the other hand, it is a general task to explore the 
structure of the electromagnetic field and the geometry 
it defines, see e.g. 

There are two main effects in the realm of ray optics 
based on the Maxwell equations: One effect is birefrin- 
gence and the other one anisotropy of the propagation of 
light Both effects are well known from the physics 
of light propagation in general media, such as in crys- 
tals, e.g.. The basics of the general formalism have been 
laid down in Q • The explicit calculations of these effects 
have been carried through to first order in these anoma- 
lous effects by Kostelecky and coworkers and by others 

H> ^ IE 1^ (for a possible birefringence caused by a 
torsion of spacetime, see 0, H, S ^3 and also ^3). In 



these approaches the first step is to confront the result 
with the possible observations of birefringence. From as- 
trophysical observations 1,], the parameters responsible 
for birefringence must be smaller than 10~^^ and, thus, 
can safely be neglected. The remaining anisotropy in 
the photon propagation is given by a symmetric second- 
rank tensor. By an appropriate coordinate transforma- 
tion, this tensor becomes proportional to the unit tensor. 
Accordingly, there is an adapted coordinate system such 
that light propagation is isotropic and defines a Rieman- 
nian metric. This is a remarkable result that may be due 
to the approximation used. In this work we show that 
this result holds exactly. That is, we show, provided we 
assume a local and linear constitutive law for the vacuum, 
that 



Maxwell equations 

- only real sols, of Fresnel eq. 

- vanish, birefringence in vac. 



pseudo- Riemann- 
ian metric . 
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II. OBSERVATIONAL AND EXPERIMENTAL 
FACTS 

As discussed, the best estimate on birefringence effects 
of the vacuum have been given by an analysis of Kost- 
elecky and Mewes Qj. Their results show that the bire- 
fringence parameter is smaller than 10~^^. This estimate 
is independent of the coordinate system chosen since it 
is an effect which cannot be transformed to zero. 

Since the resulting anisotropy can be transformed 
away, it cannot be understood as an effect solely within 
the photon sector. Thus, the coordinates used for the de- 
scription of the anisotropy experiments have to be fixed 
by some other physical process. In these experiments, 
this is realized by some solid like, e.g., the interferometer 
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arm or the optical resonator. The length of the resonator 
or of the interferometer arm is determined by the laws 
of quantum physics and of electrodynamics, see |0, 0| , 
e.g.. Consequently, the search for an anisotropy of the 
propagation of light has to be interpreted as a comparison 
between the laws of quantum physics, like the Dirac equa- 
tion, the Pauli exclusion principle, etc., and the Maxwell 
equations. The most recent experiments searching for an 
anisotropy of the velocity of light yield no effect to the 
order of Ac/c < IQ-^^ 



III. SOME PREMETRIC ELECTRODYNAMICS 

The Maxwell equations, expressed in terms of the ex- 
citations V, H and the field strengths E, B, read 



dV ^ p, dn-V = j, 
dB ^ 0, dE + B :^0. 



(1) 
(2) 



We mark the exterior derivative in 3 dimensions with an 
underline: d. The dot denotes a Lie derivative with re- 
spect to the vector field dt ■ The electric charge density 
is p, the current density j. For the formulation of the 
Maxwell equations, we use the calculus of exterior dif- 
ferential forms. We take the notation from 0, compare 
also Frankel 0, Lindell 0, or Russer ^.g.- 
The 4 dimensional form of the Maxwell equations 

dH = J , H = V-Hhdt, J^p-jAdt,{3) 
dF=0, F = B+EAdt, (4) 

shows that they are generally covariant under diffeomor- 
phisms and there is no need of a metric of spacetime 

The set of equations and Q is incomplete. What is 
missing is the constitutive law of the vacuum (the space- 
time relation). If we assume locality and linearity, then 
H = k{F), with the local and linear operator k. If we 
decompose the 2-forms H and F in their components 
(here i,j = 0,1,2,3), then H = Hij dx' A dx^ /2 and 
F = Fi-j dx' A dx^ /2. Accordingly, 



1 

2 '-y 



^ii ^ '^ii ^ kl 



with Kij'^'- 



kl 



Ik 



(5) 



Here Ky *^' is the constitutive tensor of spacetime with 
36 independent components. With the help of the con- 
travariant Levi-Civita symbol e'^™" = ±1,0, we can 
introduce the equivalent constitutive tensor density of 
spacetime. 



X 



ijkl 



1 



Ajmn kl 



(6) 



Incidentally, the covariant Levi-Civita symbol, which we 
will use below, is denoted by a circumflex: Cijmn — ±1,0. 
Since no metric is available at this stage, we have to 
distinguish these two symbols. 

Alternatively, we can express jSJ in a six component 
version, which is sometimes more convenient. In terms 



of blocks with 3-dimensional indices a, 6, 
find 



B" 



1,2,3, we 



(7) 



Obviously, 21 is the 3-dimensional permittivity matrix 
and 05 the reciprocal of the permeability matrix. The 
matrices £ and S) describe electric-magnetic cross terms 
(which vanish in Maxwell-Lorentz vacuum electrody- 
namics in Cartesian coordinates). In ((T)), for the compo- 
nents of the electromagnetic field, we took a vector-like 
notation 



V ^ V'' ib, B = B''ib, 



(8) 
(9) 



with the 3-dimensional coframe and the 2-form basis 
Ea = ef)cd ^'^ A 'd'^/2. By straightforward algebra, the con- 
stitutive 3x3 matrices 21, *B, £, S can be related to the 
4-dimensional constitutive tensor density lO by 



OaOb 
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<^bcdX 



cdOa 



■ ^acd X 



Obcd 



(10) 

(11) 

(12) 
(13) 



IV. 



QUARTIC WAVE SURFACE FOR THE 
PROPAGATION OF LIGHT 



The propagation of light in local and linear premetric 
vacuum electrodynamics is characterized by the general- 
ized Fresnel equation 



Mok^ + Mikl + M2kl + Ahko + Mi = , 



(14) 



where ko is the zeroth component of the 4-wave covector 
k. The coefficients Mi are homogeneous functions of de- 
gree i in the spatial components ka of the wave covector: 



M, ■.= M'''-''^ka,---ka^ 



(15) 



The Fresnel equation results from a solvability condition 
for a 3-vector equation W^^ki, = on the jump surfaces 
Ulii; here 

{kl 21'^'' + k^kd [e:"c + 1\ e'^'l 

+k,kfe''^^e^f''^,a) (16) 



ab 



is a 3 X 3 matrix, the determinant of which has to vanish, 
see 18]. Eq. (|15|l is valid in a special anholonomic frame 
with 'd° ^ k. 

The equation for the jump surfaces can also be ob- 
tained in an analogous way as effective partial differen- 
tial equation for the components of the radiating elec- 
tromagnetic potential after removing all gauge freedom. 
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This equation, for all initial conditions or all sources of 
sufficient regularity, should possess a unique solution in 
some future causality cone (this corresponds to a finite 
propagation velocity of the solutions). The necessary and 
sufficient condition for that is the hypcrbolicity of the 
differential operator Furthermore, the differential 

operator is hyperbolic if the corresponding polynomial 
is hyperbohc jl^. This means that ((Tl|l is required to 



possess four real solutions for feg which need not to be 
different. The condition of the hypcrbolicity or, equiva- 
lently, the condition for the existence and the uniqueness 
of the solutions, is the fundamental fact behind the par- 
ticular signature for the metric which we are going to 
derive (see also |2(|] for another example). 

The M°i - ''''s in (|15|l are cubic in the 3x3 matrices 
21, iB, £, and 2), see U 



M 

j\^abcd 



deta, 

-ibcd (21^'^ a'^"^ + 21'^'' 21"'= D^'^) , 

i 21^'^'') [{t^'df + (S/)' - {t'd + + ^/)] + (C'^c + I)c'')(2l"^"G:''<i + Drf('^2l'')") 



^de(c\ 



^ef{c^\gh\d (g 



hf 



(17) 
(18) 



(19) 
(20) 

(21) 



Computer plots of the 4th-order surface of the general- 
ized Fresnel equation H14|l have been prepared by Terty- 
chniy 21]. 

We solve ()14|l with respect to the frequency kg , keeping 
the 3-covectors ka fixed. We find the four solutions 



with 



/a - 



^O(l) 



k T 

'^0(2) 



k ^ 



k ^ 
'^0(2) 



We introduced the abbreviations 



/3 + 


7 

\/a 


-5, 


(22) 


f3 + 


7 




(23) 


/?- 


7 




(24) 


f3- 


7 


-<5. 


(25) 



a := 

/3 := 

7 := 

(5 := 



12M, 



V (b + V^) 



1 



12Mo \^ {b + ^) 
1 



- + {b + V^)' ^ 2A/2 j + , 

(26) 



4Mo 
Ml 
4Mo 



(2(^M2 - M3) - 2<5^ 



(27) 
(28) 

(29) 



a := UMqAU - 3M1M3 + , 
27 

6 := —MqMI - 36M0M2M4 

9 27 
--M1M2M3 + —MlMi 

c := 4 (6^ _ a3) . 



(30) 

(31) 
(32) 



Earlier investigations on light propagation in general 
linear media and on Fresnel-Kummer surfaces includes 
the important work of Schultz et al. and Kiehn et al. 



V. VANISHING BIREFRINGENCE 

Vanishing birefringence means that there is only one 
future and only one past directing light cone. In order 
to achieve this, one has to identi fy tw o pairs of solutions. 
There are these two possibilities 



'^O(l) ~ '^0(2) 



h ^ — h ^ 
'^O(l) ~ '^0(2) 



/3-7 = 
a = 7 - 

For the case (|33|l . the solution degenerates to 



'^O(l) ^ '^O(l) 



h ^ — h ^ 
'^0(2) ^ '^0(2) 



I.e. 
i.e. 



and for the case H34|) to 



k^ 



-\/a — S , 



0, (33) 
0. (34) 

(35) 
(36) 
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The equation 7 = 0, which is vaHd for both cases, has 
the simple solution 

This can be inserted into a and 6, but presently we don't 
need the explicit expressions. The functions a and (3 can 
be written as 



3Mf- 8M0M2 

6Mf - I6M0M2 ^ 
^ ~ 48M2 ^' 



with 



12M, 



'^0 V {b + ^) 



r + + 



(38) 
(39) 

(40) 



Since either /3 = or q = 0, we can add H38I) and H39() 
and find 



a or P 



3M2 - SAfoMa 
16M§ 



(41) 



corresponding to or to respectively. 
Hence in all cases the light cones turn out to be 



i,Ti _ I / 3-^1 - 8M0M2 Ml 



4Afn 



(42) 



Accordingly, the quartic wave surface in this case reads 

[{ko - kl){ko - k^)]^ = . (43) 
We drop the square and find 



ko 



Ml iM^ -8M0M2 



16Af2 



,, Afi /3A//2-8A//0A//2 \ ^ ,,,, 



4Afo 

Multiplication yields 



4A//n 



fco + ^V-^^\:lf°^-=0 (45) 



16Afn^ 



, 1 A//i 1 M2 1 /A//i\^ , , 

«o + o77^fco + 77#"o 77^ =0. (46) 



2 A/o 2 Afo 8 V^^o, 
If we substitute the A^^'s according to (|15|1 . we have 

g fcj/Cj :— kQ + — koka 



1 ( Ar'' M''M^\ 
8 \ ~I M^^^'-^" 



(47) 



This form is quadratic in the wave 4-covector ki and 
thus constitutes, up to a scalar factor, a Riemannian 
metric. Equation (|47|l represents our main result. It 
is clear that there is a coordinate system so that the 
metric g^^ acquires the ordinary Minkowski form: g^^ = 
diag(+l, — 1, — 1, — 1). Therefore, intrinsically it is not 
possible to have an anisotropic speed of light. 

From the condition of the existence of a unique solution 
(or from hyperbolicity) , equation 147|l has to possess two 
real solutions for any given spatial ka- As a consequence, 
the signature of the metric g^^ is (+1, — 1, — 1, — 1). Ac- 
cordingly, the signature of the metrical structure is a 
consequence of the existence of a unique solution of the 
Maxwell equations in a future causal cone for arbitrary 
sources with compact support. 

Let us look at a specific example. If we exclude, besides 
birefringence, also electric-magnetic cross terms in the 
spacetime relation O, then £ = J) = and, according 
to (UHl, M"- = 0. If we substitute this into (gZI), we find 



fco' 



2Af 



= 0. 



(48) 



It can be shown [2J| that one arrives also at this re- 
sult by only forbidding the existence of electric-magnetic 
cross terms, that is, this condition is stronger than the 
requirement of vanishing birefringence. Clearly then, for 
the Minkowskian signature we have 



M^^kah 
2M 



<0, 



(49) 



see also H47|) . The flat Minkowski spacetime of special 
relativity is a subcase of H48(l . Then, in Cartesian co- 
ordinates, Ai^"** is a constant. This is a consequence of 
the constancy of the constitutive matrices 21 '"^ and 05 6a. 
Because of O, we find = -%'"'Eb and Ha = ^baB^ . 
Thus, 



1 

Mo 



(50) 



where 1-? denotes the 3-dimensional unit matrix. If we 



substitute this into d)) to , we find M 
0, and A4""'' = (2eoVA*o) that is. 



-0! 



Af'^ 



2A4" 



(51) 



is negative, with c as the speed of light in vacuum. 

Note that the vanishing of birefringence is not equiva- 
lent to the validity of the reciprocity relation as discussed 
in 12]. 



VI. A UNIQUE LIGHT CONE IS 
INCOMPATIBLE WITH A FINSLERIAN 
GEOMETRY 

Now we would like to generalize the result obtained 
above: For all hyperbolic partial differential equations 
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a vanishing birefringence of the characteristic cones de- 
fines merely a Riemannian structure. There is no way to 
have two characteristic cones with a Finslerian structure. 
In fact, the restriction to hyperboUc partial differential 
equations is necessary for physical reasons: only for hy- 
perbolic partial differential equations one has a unique 
solution in the future half space for prescribed initial val- 
ues or prescribed source, see p^ . 

Let us now prove the above statement: Any charac- 
teristic surface is given by a polynomial of order p in the 
covector ki, which is 'normal' to the characteristic surface 



(52) 



In order to be based on a hyperbolic differential operator, 
this polynomial also has to be hyperbolic, that is, there 
should exist p real solutions = ko{ka) (see, e.g., [Tojl 



Hik) - n (A:o - k 



'O(m) 



(53) 



This specifies a splitting of the characteristic cone into p 
sheets. 

Now we want to restrict the number of cones to 2. In 
order to be able to identify an equal number of cones, 
we choose p = 2q. After the identification of the first q 
solutions and the last q solutions, respectively, we have 
as characteristic polynomial 

H{k) = g^^^^-'-^h,k,,---k,,^ 

= (fco - fco(l))'' (fco - fco(2))'' 

= [ko - {ko{i) + fco(2)) ko + fco(i)A;o(2)] , (54) 

where the two solutions fcg(x) and fco(2) are homogeneous 
functions of the spatial components ka ■ 

We differentiate this relation with respect to ka and set 
subsequently ka ~ 0. This results in fco(i 2)(fca = 0) = 0. 
For the zeroth derivative we get 



(55) 



The first derivative reads 



2ng ^ 2,1 1 j^^^ . . . ki^ 



q {kl - (/co(i) + fco(2)) fco + A:o(i)/co(2))'^ x 
f d d \ 

\dK ^ ^''^'dk'a (^0(1)^0(2)) j , (56) 

which, for fc^ ^ 0, yields 

d 



20 



(O...Oa) _ 



dka 



(fco(l) + fco(2)) 



(57) 



This can be integrated to 

fco(i) + fco(2) - -2g(0-°'')A;„ =: -2g°'^ka (58) 

(no constant must be added because the fco(m)'s are ho- 
mogenous in ka)- 



Analogously, we calculate the second derivative and 
perform the limit fc^ — > 0, 



2(2g-l)fl 



(O...Oab) 



4g(0...0a)^(0...0b) 
92 



dkadkb 



(fco(i)fco(2)) , (59) 



where we used Therefore, 



^0(1)^0(2) 



(2g- 1)0(0-0"'') -20(o-o")g("-o^) 
g^'kakb . 



(60) 



If we substitute H58(l and (|6()|l into H54|l , then merely a 
Riemannian metric shows up. 



ko - (^o(i) + ko(2)) ko + A:o(i)fco(2) 



(61) 



with = 1. No Finslerian metric does occur. The un- 
derlying metric has to be of signature ±2. Otherwise 
it would not lead, for prescribed ka , to two real solutions 
/co- Again, the metric g'^ has to possess the signature 
(+1,-1,-1,-1). ■ 



VII. DISCUSSION 

As our main result, we have shown that radiative vac- 
uum solutions of the general Maxwell equations that do 
not show birefringence define — up to a scale transfor- 
mation — a Riemannian metric. Thus, the requirement 
of vanishing birefringence automatically yields a Rieman- 
nian structure. No Finslerian metric can be introduced 
in this way. As a consequence, no intrinsic anisotropy 
in the propagation of light can be found (intrinsic in 
the sense of using merely the Maxwell equations). It 
is always possible to make a coordinate transformation 
to a locally Minkowskian frame. This applies also to 
a hypothetical higher order version of the generalized 
inhomogeneous Maxwell equation like dj{x^^'^'' Fki/2) + 
djdm{x^^"^^'' Pki) = J*. Only if non-Minkowskian coordi- 
nates are related to or fixed by other physical phenomena, 
then one may speak about an anisotropy of the speed of 
light. Such phenomena may be related to quantum mat- 
ter described by some Dirac-like equation. Accordingly, 
this anisotropy is defined only with respect to another 
physical phenomenon. 

This situation is, of course, present in current tests 
searching for an anisotropy of the propagation of light, 
like the modern tests using optical cavities ■ In these 
tests the isotropy of the velocity of light is tested with 
respect to the length of the cavity. This length is deter- 
mined by the Dirac equation but, in part, also by the 
Maxwell equations. However, it turns out that for the 
used materials the latter influence the length of the cavity 
only marginally so that the length is mainly determined 
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by the Dirac equation. Therefore, Michelson-Morley 
tests are tantamount to a comparison of the Maxwell with 
the Dirac equation. 

This result also shows that the generalized Maxwell 
equations alone cannot cover the anisotropy effects 
of light described in the kinematical framework of 
Robertson-Mansouri-Sexl . In the same way 

as in this kinematical framework, one has to make a com- 
parison between the propagation of light and a length 
standard. This length standard is given as such within 
the Robertson-Mansouri-Sexl framework. In the present 
framework of dynamical test theories, this is replaced by 
a comparison of the Maxwell and the Dirac equation. In 
this sense, one may take the framework including a gen- 
eralized Maxwell and a generalized Dirac equation as the 
dynamical replacement for the old Robertson-Mansouri- 
Sexl framework. However, one may want to go further 
to the appreciably more general standard model exten- 
sion (SME) of Kostelecky and collaborators 28], which 
contains more than a single generalized Dirac equation. 
In any case, the birefringence of light and also of Dirac 
matter waves in vacuum is truly beyond the Robertson- 
Mansouri-Sexl scheme, but is included in the SME of 
Kostelecky. 

Our main result only relies on the fourth order Fres- 
nel equation (|14|l . All propagation phenomena which 
lead to characteristic equations of fourth order lead to 
a Riemannian metric if one does not allow birefringence. 



Therefore, this also applies to the characteristics of a 
generalized Dirac equation where the 7-matrices are not 
assumed to fulfill a Clifford algebra. If the Dirac charac- 
teristics do not show birefringence, then we can conclude 
that the 7-matrices will fulfill a Clifford algebra. This 
also follows from our general result in Sec lVIl 

Furthermore, our result can also be applied to the 
WKB approximation of generalized particle field equa- 
tion as, e.g., the generalized Dirac equation p0l l29ll30| . 
As a result, one arrives at a scalar-vector-tensor the- 
ory where the dispersion relation induces a splitting of 
the mass shells according to = fcg — g°'''{pa + Oia){pb + 
ab) + o? ■ The equation of motion for the corresponding 
point particle is that of a charged particle in Riemannian 
spacetime with a position and time dependent mass. 
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